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Abstract 


This  paper  is  the  first  in  a  series  of  three,  which  analyze  an  adaptive 
approximate  approach  for  solving  (n+1)-  dimensional  boundary  value  problems 
by  replacing  them  with  systems  of  equations  in  n-dlmensional  space.  In  this 
approach  the  unknown  functions  of  (u+1)-  variables  are  projected  onto  finite 
linear  combinations  of  functions  of  just  n-variables. 

This  paper  shows  how  the  coefficients  of  these  linear  combinations  can 
be  chosen  optimally. 


Introduction 


Vi  n+T  V. 

Let  u  m  a)  x  [-h»h]  be  a  domain  in  IR  ,  and  let  u  be  the  solution 
to  some  elliptic  boundary  value  problem  on  .  We  wish  to  find  —  in  a  very 
effective  way  —  an  approximate  solution  Uapprox  satisfying  a  certain  accuracy 
requirement. 

Considering  the  special  structure  of  Q*1  ,  we  expect  that  u*  can  be 
approximated  well  by  a  linear  combination 


N 

l  ’My/h)*x.  .  y€[-h,h] 

j-0  2  2 


)N 


of  N+l  functions  {x^ }  ^  on  u> 


Methods  built  on  an  assumption  of  this 
type  and  a  projection  procedure  are  widely  used  in  engineering.  As  an  example, 
we  mention  various  theories  for  plates,  beams,  etc.  (cf.  [4,11,12]).  These 
methods  are  also  sometimes  associated  with  the  name  of  L.V.  Kantorovich  (cf.ClO]). 


Our  goal  is  to  select  the  family  of  functions  suc**  that 

i)  The  method  is  optimally  accurate  when  h  is  small  and  the  data 
sufficiently  regular. 

ii)  For  arbitrary  h  and  input  data  the  method  converges  as  N  -*>  •  . 

iii)  It  is  possible  to  derive  an  a-posterlori  estimate  for  the  error  and 
this  leads  to  an  effective  procedure  for  the  selection  of  N  . 


Another  approach,  which  has  also  been  extensively  used  in  structural  mechan¬ 
ics  and  elsewhere  to  derive  lower  dimensional  approximating  models,  is  asymptotic 
expansion  in  h  (cf.  [5,8,9]).  We  refer  to  [16]  and  references  therein  for 
various  engineering  applications  of  this  approach.  It  is  quite  obvious  that  a  method 
based  on  an  asymptotic  expansion  in  h  does  not  satisfy  the  goals  stated  above.  The 
validity  of  an  approach  of  this  type  is  dependent  on  the  smallness  of  h  ,  while  the 
actual  value  of  h  may  simply  not  be  small  enough.  But  even  for  arbitrarily 
small  h  the  approximation  can  be  unsatisfactory  because  of  rough  data.  [1] 
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contains  an  example  that  shows  that  even  for  an  extremely  thin,  simply-supported 
polygonal  plate,  the  biharmonic  equation  is  not  always  a  good  model  for  a  three- 
dimensional  plate. 

For  the  model  problem  analyzed  in  this  paper  we  restrict  ourself  to  pro¬ 
jection  in  the  energy.  We  project  the  solution  xx1  on  elements  of  the  form 


N 


l  y/h)x^ 


j-0 


03 

where  is  a  sequence  of  functions  on  [-1,1]  .  The  variable  y  ranges 

over  [-h,h]  and  the  x^’s  are  arbitrary  elements  of  some  linear  space  (e.g., 

a  function  space  on  oj  )  .  The  main  question  addressed  here  is  how  to  select 

the  sequence  ^as  alrea(*y  been  stated  that  one  of  our  goals  is  to 

obtain  optimal  rate  of  accuracy  for  small  h  .  In  theorems  3.1  and  4.1  we  prove 

00 

that  this  requirement  almost  uniquely  determines  the  sequence  {  ^ j.q  • 

In  a  forthcoming  paper  [15],  we  study  convergence  properties  of  the  method 
as  N  ■+  •  . 


The  presence  of  a  strong  singularity  in  the  data  is  reflected  in  the  fact 
that  we  must  use  a  relatively  high  number  of  functions  t|^  (i.e.,  we  must  in¬ 
crease  N  )  .  Because  such  singularities  are  often  localized,  it  seems  appro¬ 
priate  to  introduce  the  possibility  of  using  a  different  N  in  different  parts 
of  the  domain  ft*1  .  The  a-posteriori  error  estimation  and  the  problem  of  how 
to  design  an  adaptive  algorithm  that  will  produce  a  good  distribution  of  the  N's 
was  briefly  discussed  in  [14].  It  will  also  be  given  a  more  detailed  treatment 
in  a  forthcoming  paper. 
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2. 


Notation  and  the  Model  Problem 

Let  H  be  a  separable  Hilbert  space  with  Inner  product 


•  <u,u> 


1/2 


<u,v> 


and  norm 


A  denotes  a  (possibly  unbounded)  self  adjoint  linear  operator  In  H  with 
domain  of  definition  V(A)  .  Furthermore  we  assume  that  A  Is  a  strictly  positive- 
definite  operator,  i.e.,  there  exists  C  >  0  such  that 


Vu€P(A):  c|  |u|  <Au,u> 


t>(A1/2) 

<A1/2»,A1/2v> 


is  itself  a  Hilbert  space  with  the  inner  product 


<u,v>  + 


Let  M  be  a  self  adjoint  bounded  linear  operator  in  H  .  M  is  also 
assumed  to  be  a  strictly  positive- definite  operator, 

I  denotes  an  interval  on  the  real  line.  L  (I;ff)  is  then  defined  as  the 
set  of  strongly  measurable  functions  I  -*•  H  such  that  ||u(*)||  is  an  element 
of  L^(I) ,  (cf.  [6]).  The  same  goes  for  L^(I;P(A*^))  • 

We  also  need  a  Sobolev  space  of  functions  with  values  in  H  .  H^(I;H) 

2 

denotes  the  space  of  functions  I  -*■  H  such  that  u(*)€L  (I ;H)  and 
d  2 

^u(*)€L  (I;H)  ,  (df  .  [2]).  The  derivative  here  is  taken  in  the  distributional 
sense.  H^(I)  denotes  the  standard  Sobolev  space  on  I  . 

Assume  a  and  b  are  real  valued  functions  in  L^Q-l,!])  such  that 


ao  Vy€[-l,l]  , 

b0  £  b(y)  Vy  € I-i.U  , 

for  some  constants  a  >  0  ,  b  >  0  . 

o  o 
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a^  and  b^  CL^d-h,!*])  are  then  defined  as 

ajj(y)  -  a(y/h)  Vy€[-h,h]  , 

bh(y)  ■  b(y/h)  Vy€  t-h,h]  . 

By  ("J~)  we  denote  the  differential  operator  -  ^  (a^ 

Let  f  and  g  be  two  arbitrary  vectors  from  H  .  We  consider  the  fol¬ 
lowing  model  problem 


(1) 


ph<w)M"h  +  bh Auh 

.  d  u  h 

*h  d7Mu 

.  d  V  h 

\  d7Mu 


g 

f 


in  ]-h,h[  , 

for  y  =  h  , 
for  y  =  -h 


(Other  boundary  conditions,  e.g.,  Dirichlet  conditions,  could  just  as  well 
have  been  chosen;  we  could  also  consider  the  inhomogeneous  problem.  The  above 
selection  was  simply  made  for  convenience.) 

Before  we  proceed  any  further,  let  us  give  a  simple  example. 


Example 

Let  u  be  a  domain  in  lRn  with  a  Lipschitz  boundary.  As  H  we  take 

2 

L  (u>)  .  Let  A  be  the  Friedrichs  extension  (cf.  [13])  of  the  operator 

Ol 

-div  c(x)  grad  defined  on  a  subspace  of  H  (w)  .  (x  denotes  coordinates  in  oj) . 
c  is  a  function  in  L*(u>)  such  that  3c  >0  with  c  <  c(x)  w  x€(»)  . 

If  we  take  a  •  b  and  let  M  be  the  operator  of  multiplication  by  c(x) 


the  problem  (1)  becomes 


div  dh(x,y)  grad  u  *  0  in  cu  x  ]-h,h[  , 


^(x.y)  uh  *  g(x)  for  y  -  h  , 
^(x.y)  ^  “h  "  *00  for  y  -  -h 


u  ■  0  on  3(j  x  [-h,h] 

Here  *  b(y/h)*c(x)  ,  and  div  and  grad  are  taken  with  respect 

to  the  n+1  coordinates  (x,y)  . 

The  precise  formulation  of  (1)  is 


h  —  ,,1 


uw€HJ‘([*-h,h];W)  OL*([-h,hJ;  Z?(A1/2))  , 


V»  ,v)  -  <g,v(h)>  -  <f,v(-h)> 


Vv€H1([-h,h];H)ni.2([-h,h];  P(A1/2)) 


where  8^  denotes  the  bilinear  form 


h  n 

8h(u,v)  -  j  ah  <M1/2  ^  u,  M172  ^  v>  dy  +  J  bh  <A1/2u,  a1/2v>  dy 


If  H^([-h,h];H)nL2([-h,h];P(A^72))  is  endowed  with  the  natural  norm 


n  n  . 

"  tj  I  u(y)|  |2dy  +  |  |  |Al72u(y)  |  |2dy] 


then  It  is  not  difficult  to  prove 


Proposition  2.1 


There  exist  C^t  C£  (independent  of  h)  such  that 

(i)  |8h(u,v)|  <  Cjllulll-IIMU  , 

(**■)  I  I M I  I2  1  c2  |8h(u,u) | 

Vu,v  CH^([-h,h];W)  riL2([-h,hJ;P<A1/2)) 

Also  one  has 

Proposition  2.2 

If  x€H  and  yo€t-h,h],  then 


A:  v  ■*>  <x,v(yQ)> 


is  a  continuous  linear  functional  on  Hx([-h,h];H) 


Proof 


From  the  definition  of  H^([-h,h];H)  it  follows  that  if  v(*)€  H^([-h,h];  H) 


and  X  is  a  continuous  linear  functional  on  H  ,  then 


X(v(-))€Hx([-h,h])  and 


n  n 

[j  |X(v(y))|2dy  +  j  X(v(y))|2dy] 


a  n  j 

i  IUII'  [  j  I  lv(y)  1 1 2dy  +  |  11^- v(y)||2dy] 
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(| | • | | '  is  the  norm  in  H’)  . 

Hence,  with  A(*)  *  <x,  •>  we  get 


<x,v(* )>  €H*([-h,h])  and 

h  h  1/2 

[J  |<x,v(y)>|2dy  +  |  1^  <x,v(y)>  1 2dy] 

-h  -h 


< 


|v(y)  |  |2dy  + 


a  ,  11/2 

J  1 1^  v(y)l  l2dy] 

-h 


Using  this  last  estimate  together  with  the  standard  trace  theorem,  we 


finally  get  that  A  is  a  continuous  linear  functional  on  H  (f-h,h];H) 


□ 


Propositions  2.1  and  2.2  immediately  give 
Proposition  2.3 

The  problem  (2)  has  a  unique  solution. 


3.  The  Direct  Result 

We  first  define  exactly  what  is  meant  by  a  dimensionally  reduced  solution 
00  2_ 

to  (2)  .  Let  }j_qCH  ([-1,1])  be  a  given  sequence  of  linearly  independent 

functions. 

Definition 

The  dimensionally  reduced  solution  u^  of  order  N  is  the  projection 
of  u1  onto  the  space 


-  {  l  *  (y/h)x.|x.€P(A1/2)  ,  j  “  0, . . . 
j-0  3  J  -> 


The  projection  is  with  respect  to  the  inner  product  8.  (u,v)  . 

n 

The  main  result  of  this  section  is  the  next  theorem,  which  suggests  a  way 
of  choosing  a^so  S*-ves  an  estimate  of  |  |  |u^-u^|  |  |  in  terms  of 

powers  of  h  . 

Let  F  denote  the  differential  operator 


k"1  d  J_ 

dy  dy 

In  the  notation  of  the  previous  section  P  *  -b  ^  P^(x^)  . .  P  is  considered  as 
an  operator. 


L2([-l,l]) 2  P(P)  ♦  L2( [-1,1] ) 


N(P^)  denotes  the  null  space  of  the  operator  P*  ,  0  £  i  .  P°  *  I  (identity), 


It  is  easy  to  see  that  N(P^)£  H2([-l,l])  for  all  i  . 
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Theorem  3.1 

oo 

There  exists  a  sequence  of  linearly  independent  functions  ,  with 


N( P1)  -  span  1  1  1  » 

that  has  the  following  property: 

For  any  integer  N  >_  0  and  for  any  given  set  of  vectors 
-1  N 

£f  g€V((AM  )  )  there  exists  a  constant  (independent 
of  h)  such  that 


,  h  h 

lu  ~U2N 


2N+1/2 


Remarks 

00 

The  sequence  depends  on^Y  on  the  operator  P 


It  is  also  clear 


that  ^ j*o  *s  not  uni*uely  determined  by  Theorem  3.1.  Any  other  sequence 
,  with  span  *  span  ^ j„o  ^ *  9  could  have  been  used. 


The  rest  of  this  section  is  devoted  to  proving  Theorem  3.1.  To  do  so  we 
need  a  couple  of  auxiliary  results. 

By  changing  variables  to  [-1,1]  and  introducing  uh(y)  *  u*(y*h)  for 

Mi 

“1  i  y  i  1  we  transform  (2)  into  the  following  equation  for  u  : 


(3) 


uh€H1<C-i.i3;tf)  nL2([-i,i];P(A1/2»  , 

8h(uh,v)  -  <g,v(l)>  -  <f,v(-l)>  , 

Vv€H1([-l,l];H)  DL2(C-1,1];P(A1/2)) 
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Here  the  bilinear  form  8^  is  given  by 

^(u,v)  -  h"1  |  a<  M1/2  H1^2  dy  +  h  J  b  <A1/2u,  A1/2v>  dy  . 


-1 


-1 


Now  let  us  define  the  sequence  by  the  following 

equations: 

For  any  v€H1(C-l,l]) 


(4°) 


fa^dvd 
J  dy  dy  dy 
-1 


0  , 


(41) 


1  A,° 

dtp 


f  a*l  dv  f  K)o 

J  a  W  &y  dy  J  b*o 

-1  -1 


|V  dy  =  v(l)  , 


and  for  j  ^  2; 


(4j) 


1  % .  o 

dtp 


j a  %  *•> +  J  S-i 

-1  -1 


v  dy  *  0 


1  00  1 

The  sequence  ([-1,1])  is  defined  by  the  same  system  of  equa- 
tions,  the  only  difference  being  that  in  the  right  hand  side  of  (4^)  v(l)  is 
replaced  by  v(-l)  . 


Lemma  3.1 


Let  j  denote  an  integer  >  0  . 


i  t  j 

The  equations  (4J)  0  j  <_  Jq  determine  the  set  un^<Jue^y 


up  to  a  constant  in  ip.  >  -£*0,1. 

^o 


r 


S: -%  '-s 


Proof 

It  is  sufficient  to  prove  the  statement  for  t  ■  o  . 

The  case  j  •  0  is  obvious.  We  proceed  by  induction.  Hence  assume 
°  j  -l 

that  the  equations  (4^)  0  _<  j  <_  jQ  -  1  determine  uniquely  modulo 

a  constant  in  . 

1  "1  • 

Jo 

Consider  the  equations  (4^)  0  j  j  .  According  to  the  induction 

o  j  -1 

hypothesis,  q  is  determined  uniquely  up  to  a  constant  in  ij>  ,  . 

3  j  1  o 

Choosing  v  »  1  in  equation  (4  °)  we  derive  the  value  of  /  b  ^  dy  , 

•1  j 

which  means  that  ij/°  is  completely  determined.  The  equation  (4  °)  is 

3o  d  d  1  o 

now  nothing  but  a  Neuman  problem  for  -  ^  a  ^  ,  and  since  /  b  ^  v  dy 

,/JOv 

is  equal  to  the  right  hand  side  of  equation  )9  for  any  constant  v  ,  this 
has  a  solution  that  is  unique  modulo  a  constant.  This  proves  that  is 

determined  uniquely  by  up  to  a  constant. 


q  oo  ^  00 

Because  of  the  way  the  two  sequences  an<*  j“0  are  constructe^ 

we  also  have 


Lemma  3.2 


For  any  1  0 


NCP1)! 


<£span{t|yi|>j  N(Pi+1) 


Proof 


The  lemma  Is  clearly  true,  if  for  any  i  >_  0  we  can  prove  the  more  detailed 


statement: 
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A/(Pi)^span{^,^>jj.0^W(P1+1)  , 

o  1 

and  there  exists  a  linear  combination  ♦  ”  °^i+l  +  ®^l+l  suc^ 
that 

spanf'l'j  f  }j_q  &  span{i|<>  *  N(Pi+1) 

The  validity  of  the  above  statement  is  easily  checked  for  1*0,  and  we 
proceed  by  induction.  That  is,  we  assume  the  statement  to  be  true  for  i*k>0 
From  the  way  *^+1  an<^  ^k+1  ar®  constructed  it  follows  that 

*kn-  *k«*M<plt+2,V'<pW’1)  • 

and  because  of  the  induction  hypothesis  we  then  get 

N(Pk+1)^span{i|»°,t|»j}^  . 

Now,  adding  the  two  vectors  *£+1  *nd  ^  ,  we  cannot  increase  the  dinen- 
sion  by  more  than  two.  The  fact  that  the  codimension  of  N(Pk+1)  in  N(P  ) 
is  2  ,  together  with  the  strict  inclusion 

span{i|»°,0j}k-o^N(Pk+1)  , 

hence  shows  that 

span  { t|»° ,  W  (pk+2  ) 
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That  means  we  have  proven  the  first  part  of  the  extended  statement  for 
i  -  k+1  . 

Concerning  the  second  part  we  consider  the  linear  combination 

“4+2  +  6  4+2  * 

From  the  construction  we  know  that 

p(a4+2 +  g4+2> "  a4+i +  8 4+i  • 

Together  with  the  induction  hypothesis  this  tells  that 


+  «*i+2«N<>'fe+2)- 


On  the  other  hand,  “^+2  +  ^^1+2  cannot  an  element  of  spant^j 
If  so,  we  would  have  by  application  of  the  operator  P  that 


a4+i  +  84+i€sPan{4’4}d-o  ’ 

which  contradicts  the  induction  hypothesis. 

This  finally  proves  that  “^^.2  +  8  4+2  to8ether  with  sPan^j*4^^0 
spans  all  of  N(pk*^)  . 

Thus  the  second  part  has  been  established  for  i  ■  k+1  . 

□ 

In  the  proof  of  Theorem  3.1  we  shall  also  use  the  following  density  result. 


jrtfcSfBtv-.  ..  —c-'-j, 
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Lemma  3.3 

The  set  V  defined  as 

J 

V«{£  v.  x.  |  J63N,  x.  £P(A^2)  and  v.  €H*([-1,1])  for  0  <_  J  <  J} 
j-0  J  J  J 

Is  dense  in  1];H)  D  L2([-1,1];P(A1/2)  . 


Proof 

The  lemma  is  clearly  proven  If  one  shows 
(i)  V  is  dense  in  H^C-l.l]  ;I7(A1/2))  ; 

(ii)  H1([-1,1];P(A1/2))  is  dense  in  H1(C-l,l] ;H) nL2([-l,l];P(A1/2))  . 

Let  us  start  with  (i).  Asstime  u€VA  (A  in  H*(  [-1,1  ];P(A*^2))  )  ,  that 


<A1/2  ^ 


du 

dy’ 


A1/2„> 


-1 


dvd  J 

-7-  dy  + 

dy  3  J 
-1 


<A1/2», 


p}l2-x> 


v  dy 


VvfH^-l,!]),  x€0(A1/2) 


Since  A  is  closed,  this  yields 


J  dy 
-1 


<A1/2u,  A1/2x> 


%iy  +  \ 

-1 


<A^2u,  A^2x>  v  dy  ■  0 


Vv  €H^([-l,l]) ,  x€P(A^2),  and  hence 


<A1/2u,  A^2x>  -  0  V*€0(A1/2)  , 


■&£&  f 
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or 

<u,Ax>  *  0  Vx€P(A) 

Because  of  the  fact  that  A  is  invertible,  we  know  that  A(P(A) )  -  H  .  From 
the  previous  equation  we  thus  get 

u  *  0 


In  summary,  we  have  now  proven  that  V1  ■  0  .  This  immediately  yields 


(i)  . 

1/2 

To  prove  (ii),  let  dE^  be  the  spectral  measure  of  A  (cf.  [7]). 
If  u€H^([-l,l];H)  nL^([l,l];P(A^^))  ,  then  u  can  be  written  as 


for  k  -*■  «  ,  because  of  Lebesque's  monotone  convergence  theorem.  In  the 
same  fashion 

|  I |A1/2(u-uk) | |2dy  - 

-1 

1 

|  |  X2d| |Exu||2dy  -  0 

-1  a(A1/2)  n]k,«[ 

for  k  «♦*  «  .  That  is,  we  have  proven 

1 1  lu"uiJ  1 1  -*■  °  f°r  k  **■  °°  t 
which  immediately  gives  (ii) . 

□ 

Finally  we  are  now  ready  for  the 

Proof  of  Theorem  3.1 
00 

Choose  such  that 

spanUj}2^1  -  N(PA)  Vi  >  1  , 

and  span{ipj}2f0  -  span{i|>°,t|>j}j-0  Vi  >  0  • 

This  is  possible  because  of  Leona  3.2.  The  ij^'s  chosen  this  way  obviously 
have  the  first  property  stated  in  Theorem  3.1  • 
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For  any  N  0  and  any  pair  f,g€P((AM  )  let 

SN  “  l  h"1+2j(i|»°(y/h)M’1(AM'1)j'1g  -  ^(y/h)M"1(AM"1)j~1f ) 
j-0  3  3 

U  U 

It  is  clear  that  .  Because  of  Proposition  2.1,  it  follows  that 

there  exists  a  constant  C  (independent  of  N  and  h  )  such  that 


uh-u$1)|  1 1  i  0  Inf  II  |„"-v|  1 1  <  C|  I  |«"-sj|  1 1 


h  „h, 


v  €V 


2N 


Ue  now  estimate 


h  eh 
°  -s» 


By  a  change  of  variables  and  introduction  of  the  bilinear  form  8^  it 
follows  that 


with  sjl  -  l  h'1+23(t|»?(y)M_1(AM_1)3‘1g  -  ^(y)M_1(AM_1)3_1f )  (and  C  independent 
N  j-0  3  3 

of  N  and  h  ) . 

Let  us  now  consider 


«  Ah  ^h  . 

8h(u  -SN,v) 


where  v  has  the  form  w(y)#x  ,  x€V( A^2)  and  w€H^([-l,l])  .  We  know  that 


S^-sJ},*)  -  w(l)  <g,x>  -  w(-l)  <f,x>  -  fi^(sjj.v) 


Concerning  the  last  term  we  have 


-  ^^h“1+2:J8h(ij»®M"1(AM"1);J"1g,v) 


-  I^h”1+2^8^(i|»^M-1(AM_1)^”if  ,v) 


o,® 


Now  because  of  the  properties  of  the  sequence  {iK> 


j  j*0 


I  h‘ 1+2 j  B  (*  °M~ 1 ( AM~ 1 ) j  ~ 1 g , v) 

j-0  “  3 


1  ...o 


N  | 

+  l  h2^  btp°wdy  <  (AM_1)^g,x  > 

j-o  ix  3 


»«>  <8,*>  -  h2»  j  a  fjSti  £  4,r  <  (AH-1)",.,  >  . 


dy  dy 


-1 


Similarly 


I  h'1+2j8  (i|<}ll'1(AM’1);J'1f,v)  - 
j-0  h  J 


w(-l)  <f,x>  -  h 


1  dtp,1 


2N  r  a*Wl  dw 

J  dy  dy 
-1 


dy  <  (AM_1)Nf,x  >  , 


i  Jh  v 
8h(u  ‘SN*V) 


2N 


-1 


*  <  IT  ’  %  >  dy 


so  that  altogether 


•  *9  ' 


for  any  v  of  the  form  w(y)*x  ,  x€t?(A1^2)  and  w€H1([-l,ll)  .  rN  is 

given  by 


rH  ■  Cl<*M~1>’«  -  Cl(AM'1)Nf  • 


Now,  because  of  continuity  and  Lemma  3.3,  we  get  that 


«  Ah  £h  . 
8h(u  -SN,v) 


,  2N  f  „  drN  dv  ^  . 

h  J  a  "  17’  d7  >  dy 


Vv€H1([-l,l];«)nL2a-l,l];P(A1/2))  . 


Using  Schwarz's  inequality,  we  see  that  the  right  hand  side  is  bounded  by 


CN  h2N+^'//2[8^(v,v)  J^^2  (CN  independent  of  h). 


As  a  consequence  of  this  it  follows  that 


Ah  £h  Ml  X.l/2  „  .  2N+1/2 

[VU**SN’  U"SN)]  —  CN  h 


and  hence  |  |  |u^-sjj|  1 1  £  h2N+^^2. 


4. 


The  Inverse  Result 


In  Theorem  3.1  we  examined  a  particular  choice  of  the  functions  ifrj 

Q%T  Li  /O 

The  approximation  error  was  of  order  h  using  the  2N+1  functions 

r.  -.2N 

lVj-o  • 

The  goal  of  this  section  is  to  prove  that  the  previous  choice  of  functions 

00 

was  by  no  means  arbitrary.  That  sequence,  or  any  other  sequence 
2i  2i 

with  span  ■  span  ^or  every  *  »  c^e  on^y  one  that  gives 

this  order  of  approximation. 

We  formulate  this  as 


Theorem  4.1 

Let  N  and  K  be  two  non-negative  integers  and  f  and  g  two  linearly 

00 

independent  elements  of  H  .  Let  se<luence  introduced  in 

Theorem  3.1  and  u*1  the  solution  to  (2)  of  section  2. 

If  {*  a  set  ®le®ents  of  H^([-l,l])  with  the  property  that 


inf  |||uh-v|||  -  «(h“xt2H-3/2>-1/2,) 


v  €w£ 


where  wj^  denotes  the  set 


K 


{  l  *,(y/h)x,  I  €V(A112)  0  <  j  <  K  } 

j-0  2  2  2 


then 


2N  K 

span{^}j_Q  C  span{<|»j}j_0 


Theorem  4,1  is  actually  a  little  stronger  than  just  an  inverse  of  Theorem 

3.1.  Let  N  be  >  1  .  Theorem  3.1  then  says  that  with  the  2N+1  functions 

2N  2N+1/2 

we  can  obtain  an  error  of  order  h  .  But  Theorem  4.1  tells  us 

2N™3 / 2 

that  even  if  we  are  satisfied  with  an  error  of  order  o(ti  )  ,  we  still 

2N 

have  to  use  all  the  functions  „ 

J  j-o 

Proof 

—1  N 

We  can,  without  loss  of  generality,  assume  that  f,  g£0((AM  )  ) 

Otherwise  we  replace  u1 9  f  and  g  by  u*1  *  (HA  ^)Nu^ ,  f*  «  (MA  ^)^f  and 

g*  ■  (MArl)Ng,  which  obviously  satisfy  the  assumptions  of  the  theorem. 

Define  sjj  ,  $>^  as  in  the  proof  of  Theorem  3.1,  i.e., 

N 

sj  -  l  h~1+2j  -  *}(y/h)M"1(Alf1)J“1f) 

j-0  2  2 

and  *  sjj(yh)  . 

Then  because  of  Theorem  3.1  and  the  assumption  of  this  theorem  we  see  that 
^  v^€  such  that 

|||s|j  -  vh|||  -  a(h««ax{2N-3/2,  -l/2}) 

By  a  change  of  variables,  this  yields 
1 

«)  |  11^-  Jh)||2d,  ■  0)>  , 

-1 

1 

(ii)  J  I  |a1/2(^  -  vh)  I  |2dy  -  0(hBax{4N-4’-2}> 

>1 

^h 

The  function  v  has  the  form 
£ 

I  ♦^(y)vi  ,  with  vicPCA1^2)  for  0  <  j  <  K  . 

J-0  22  2 


WBHHMI 


1 
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Now,  multiplying  (ii)  through  by  h2  ,  we  get  because  of  the  form  of 

N 


and  v  that 


|  |^(y)A_1/2g  -  ^(y)A_1/2f  -  h  £  *  (y)A1/2v* |  1 2dy  -  0 

1-0  J  J 


K 


1/2 hM2. 


-1 


for  h  -*■  0  .  Since  for  a  fixed  K  the  set 


{  l  y)x  I  X  €H} 


j-o 


J'"  J  1  j 


is  closed  in  L  ([-1,1];  H)  ,  we  get 

,o.  ,  -l/2„  ,/.l/..\ *“l/2< 


K 


4»0(y)A“  g  -  ^Q<y)A  '  f  -  l  <K(y)x 

j-0  J  J 


The  fact  that  f  and  g  are  linearly  independent  implies  that  so  are 


.-1/2.  ,  -1/2 
A  f  and  A  g  .  Hence 


span  {^q,  <I»J}  C  span 


This  proves  Theorem  4.1  for  N  -  0  .  If  N  >_  1  we  are  not  yet  finished. 
In  this  case  we  proceed  by  induction,  i.e.,  we  assume  it  has  been  proven  that 

span  {<|»°,  C  span  Uj>j.0 

for  0  <  m  <  N-l 


Rearranging  (i)  and  dividing  through  by  h^®*2  we  get 
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-  h"2*"1  l  vj|  |2dy  -  0 

J»0  7  J 


for  h  -*■  0  . 


(Here  we  use  that  4®4-2  <_  4N-2) . 

Because  of  the  induction  hypothesis  we  know  that 


and  now,  using  that  for  a  fixed  K  the  set 


K  d* 

{jL  d7  '  XJ€H> 


is  closed  In  ,  we  conclude  that 


d*i»H  M~1  AAM-lvin  d^nH-l  -1,  -l.mf  r  d*1 

”3F“  (  )  8  ‘  ~W~  M  (  }  f  "  I  IT  xj  • 

The  fact  that  f  and  g  are  linearly  independent  implies  that  so  are 
M~1(AM~1)mg  and  M_1(AM_1)mf  .  Hence 

dCi  d*in ,  ,f^.K 

~w~  ’  “dy  '  Ss’,an{ir)j-o  • 

From  the  way  ^  and  <f>Q  are  constructed  it  is  easily  seen  that 


span{t|»Q,^Q}  ■  {constant  functions}  , 


and  this  together  with  the  induction  hypothesis  allows  us  to  conclude  that 


span{  ,  *  j  }£j£span{fj  >*.Q 
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The  induction  proof  is  now  finished,  and  finally  we  get  because  of  the 

2^j 

definition  of  • 

span{i|>j  }j^0  *  span{i|»°  ,i|)j  CspanU^  }j„q 

□ 

By  a  slight  variation  of  the  preceding  arguments  we  could  prove  the  fol¬ 
lowing  version  of  Theorem  4.1,  for  the  case  where  f  and  g  are  not  linearly 
independent,  e.g.,  f  ”  <*'g 


Theorem  4.2 

Let  N  and  K  be  non-negative  integers.  Let  f  *  erg  >  with  g€H\{0}. 
Let  ^  be  t^e  se<luence  introduced  in  Theorem  3.1. 

If  {*  }*  0  is  a  set  of  elements  of  Hl«-l,l])  with  the  Property  that 


inf,  |  |  |  u  -v 


0(h2N-!/2) 


v€W 


K 


then 


IN  K 

span{t|»j  -  cupj  >j.0  —  sPan{<(,j  }j-o 
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5.  Numerical  Examples 

Consider  the  problem 


div(a  grad  un)  -  0  in  ]0,l[x]-h,h[  , 

a  uh  -  g(x)  for  y  »  h  , 

a  j?  uh  =  -g(x)  for  y  •  -h  , 


u  *  0  for  x  *  0  and  x  *  1  , 


,  x  fa+ 

for 

y  >  0 

with 

\a. 

T 

for 

y  <  0 

• 

(a ,  and  a1  are  two  positive  constants) .  This  problem  clearly  falls  within 

T  • 

the  framework  of  our  model  problem.  Simply  choose 


b(y)  •  a(y) 
c(x)  •  1  and 


ft 


for  y  >  0 
for  y  <  0 


A  -  -(jJ  with  P(A)  -  ^([O.lDn^ao.l]) 


and  H  -  l/([0,l])  • 

The  operator  P  In  this  case  Is  given  by 


«  -Id  d 
P  •  a  -t—  a  -r— 
dy  dy 


Define  Uj  >T.0  SL  as  follows 


♦n  “  1  *  My)  •  a_1y  » 


♦2j(y)  •  |  for  J  >  1 

-1 
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and 

y, 

<i>2j+i(y)  ■  a  1  I  ^2j^t^dt  for  3  1 1  • 

-1 

Here  denotes  the  Legendre  polynomial  of  degree  k  . 

It  is  not  difficult  to  see  that  with  this  definition 

span{fj}j=0  »  spantyj}^  Vi  >.  0  , 

00 

where  is  the  sequence  introduced  in  Theorem  3.1. 

As  before  ujj  denotes  the  dimensionally  reduced  solution  of  order  N  . 
u^  has  the  form 

h  N 

*  I  <f>.(y/h)u  (x)  , 

j-0  3  J 

01  M 

with  Uj  €H  ([0,1])  .  The  vector  U  *  ^uj^j*0  ds  t*'e  9°iuti°n  to  a  two-point 
elliptic  boundary  value  problem 

,  2 

-h  K  (~)  U  +  1/h  L  U  »  F 

U(0)  »  U(l)  -  0  . 

The  matrix  L  is  diagonal,  and  the  matrix  K  has  a  band  structure.  Both  L 
and  K  are  Independent  of  h  .  We  solve  this  problem  numerically  by  expanding 
U  in  its  Fourier  series,  only  maintaining  a  finite  number  of  terms.  Since  we 
are  interested  in  studying  the  error  Introduced  by  the  dimensional  reduction, 
we  maintain  a  very  high  number  of  terms.  The  graphs  shown  here  were  computed 
using  400  Fourier  coefficients.  This  ensures  that  the  error  introduced  by 


discretization  can  be  neglected  compared  to  the  error  introduced  by  dimen¬ 
sional  reduction. 


Let  us  start  with  the  case  a+  *  a^  -  1  and  g(x)  *  tt/4  .  Figures  1,  2 
and  3  show  the  energy  error  as  a  function  of  h  by  dimensional  reduction  of 

cp 

order  0,  2  and  4  respectively.  (Note  that  the  functions  ^orm  a 

basis  for  the  polynomials  in  this  case.)  Using  interpolation  by  the  K-method 
(cf.[3])  we  know  that  g(x)  €  (H,V(A))  a#.  An  application  of  Theorem  3.1 

hence  gives  the  following  conclusions. 

1)  The  energy  error  is  of  order  h  with  dimensional  reduction 
based  on  . 

2 

li)  The  energy  error  is  of  order  h  with  dimensional  reduction 
based  on  <frg,  and  $2  • 

2 

iii)  The  energy  error  is  of  order  h  with  dimensional  reduction 
based  on  {^,  0  £  j  ^  4} 

Figures  1,  2  and  3  illustrate  the  sharpness  of  the  theoretical  results. 

Comparing  figures  2  and  3  we  see  that 

The  requirement  f ,g €0((AM~*)N)  in  Theorem  3.1  is  essential. 

If  f  and  g  are  not  sufficiently  smooth,  higher  order 
dimensionally  reduced  models  will  not  improve  the  asymptotic 
order  of  approximation  as  h  0  . 

Now  consider  the  case  where  a+  +  a+  .  In  our  computations  a+  ■  1  , 
a+  -  2  and  again  g(x)  •  ir/4  .  Figure  4  compares  two  different  dimensionally 
reduced  solutions.  For  one  of  the  dimensionally  reduced  solutions  the  poly¬ 
nomials  of  degree  <_  2  have  been  chosen  as  basis  functions  in  the  y-dlrectlon. 

For  the  other  dimensionally  reduced  solution  the  "special"  functions  4q» 

A 


Fig.  4:  Energy  error  xlO  as  a  function  of  h  , 

—  • —  • —  using  polynomials  of  degree  <_  2, 

-X - X-  using  the  "special"  functions  ♦,  and  4 
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and  $2  »  introduced  earlier  in  this  section,  have  been  used.  (Note  that  the 
"special"  functions  are  piecewise  polynomials  in  this  case.) 

Applying  Theorem  3.1  we  expect  that 

2 

The  energy  error  will  be  of  order  h  with  dimensional 
reduction  based  on  the  functions  <|>g,  ^  and  • 

Since  the  "special"  function  is  the  constant  *  1  ,  which  of  course  is  a 

polynomial  of  degree  <_  2  ,  we  also  expect  that 

The  energy  error  will  be  of  order  h  with  dimensional  reduction 
based  on  the  polynomials  of  degree  <_  2  . 

From  Fig.  4  it  is  again  evident  that  there  is  a  very  good  agreement  between 
the  theory  and  the  computational  results.  Specifically  it  is  seen,  by  comparing 
figures  2  and  4,  that 

If  the  dimensional  reduction  is  based  on  the  "special"  functions 
00 

*  t^ien  t*'e  asymptotic  behaviour  of  the  energy  error  is 
independent  of  the  regularity  of  the  solution  u*1  across  the  line 
y  ■  0  . 

A  feature  very  similar  to  this  is  well  known  for  optimally  constructed  Finite 
Element  meshes. 
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